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Abstract 

When two non-interacting plane impulsive gravitational waves un- 
dergo a head-on collision, the vacuum interaction region between the 
waves after the collision contains backscattered gravitational radiation 
from both waves. The two systems of backscattered waves have each 
got a family of rays (null geodesies) associated with them. We demon- 
strate that if it is assumed that a parameter exists along each of these 
families of rays such that the modulus of the complex shear of each 
is equal then Einstein's vacuum field equations, with the appropriate 
boundary conditions, can be integrated systematically to reveal the 
well-known solutions in the interaction region. In so doing the mys- 
tery behind the origin of such solutions is removed. With the use of 
the field equations it is suggested that the assumption leading to their 
integration may be interpreted physically as implying that the energy 
densities of the two backscattered radiation fields are equal. With the 
use of different boundary conditions this approach can lead to new 
collision solutions. 
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1 Introduction 



This paper is a study of the space-time describing the vacuum gravitational 
field left behind after the head-on collision of two plane impulsive gravita- 
tional waves. The known exact solutions are the classical solution of Khan 
and Penrose [0 and its generalisation by Nutku and Halil 0. No details 
of the derivation were given by Khan and Penrose. The Nutku and Halil 
solution was obtained using a harmonic mapping technique. The latter solu- 
tion was subsequently rederived by Chandrasekhar and Ferrari || (to quote 
from ||: "This paper is addressed, principally, to a more standard deriva- 
tion of the Nutku-Halil solution than the one sketched by the authors.") 
who developed an Ernst-type formulation for vacuum space-times admit- 
ting two space-like Killing vectors. They demonstrated that "in some sense, 
the Nutku-Halil solution occupies the same place in space-times with two 
space-like Killing vectors as the Kerr solution does in space-times with one 
time-like and one space-like Killing vector" . Thus the origin of the solution 
is still quite mysterious. Of course the task of solving Einstein's vacuum field 
equations with appropriate boundary conditions for the space-time after the 
waves collide is mathematically very complex. On the other hand the physi- 
cal picture, at least up to the appearance of a curvature singularity, is quite 
simple: two non-interacting plane impulsive waves undergo a head-on colli- 
sion and the interaction region afterwards contains backscattered radiation 
from both waves, neither of which remain plane. Our aim in this paper is 
to introduce a simple assumption based on this picture which provides a key 
to the origin of the exact solution describing the collision of two completely 
arbitrary plane impulsive gravitational waves. 

The backscattered radiation in the interaction region of space-time be- 
tween the histories of the waves after the collision determines two intersecting 
congruences of null geodesies. These are the 'rays' associated with the two 
systems of backscattered radiation. Both congruences have expansion and 
shear. The ratio of the expansions of each congruence is immediately deter- 
mined from Einstein's vacuum field equations and the boundary conditions. 
The shear of each congruence (the modulus of a complex variable in each 
case) depends in general in a simple way on the choice of parameter along 
the null geodesies (in the sense that a change of parameter induces a rescal- 
ing of the shear). We assume that a parameter exists along each of the two 
families of null geodesies such that the shear (i.e. the modulus of the 'complex 
shear') of each congruence is equal. This is the only assumption made apart 
from the usual assumption of analyticity of the solution of Einstein's field 
equations |EJ. We show that it implies, together with the field equations, 
an equation that could be interpreted physically as saying that the energy 
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density of the backscattered radiation from each wave after collision is the 
same. In addition we demonstrate how this assumption leads to the complete 
integration of the vacuum field equations. 

The outline of the paper is as follows: In section 2 the collision problem is 
formulated as a boundary-value problem. At this stage, to make the paper 
as self-contained as possible, reference is made to an Appendix A giving a 
brief summary of the construction of a plane impulsive gravitational wave 
in the manner of Penrose H. The backscattered radiation fields, existing 
after the collision, are introduced in section 3. The assumption central to 
this study is also introduced in this section and a physical implication is 
explored. Although the key assumption is simple the full integration of Ein- 
stein's vacuum field equations emerging from it is still quite complicated and 
this is described in section 4. The complications arise because the incoming 
waves are not in general linearly polarised. For readers who do not want 
to work through section 4 the considerably simpler case of linearly polarised 
incoming waves is treated in Appendix B. It is shown there that our basic 
assumption leads to the Khan and Penrose 0] solution. Finally in section 5 
the results of section 4 are summarised and contact is made with the known 
exact solutions 0,0]. 

With the use of different boundary conditions to those employed here 
the approach of this paper can lead to new collision solutions of the field 
equations. The present authors have published a new collision solution of 
the Einstein-Maxwell field equations in []5J derived using the ideas described 
in the present paper. 



2 The Boundary— Value Problem 

The line-element of the space-time describing the vacuum gravitational field 
of a single plane impulsive gravitational wave having the maximum two de- 
grees of freedom of polarisation may be written in the form (see Appendix 
A) 

ds 2 = -2 d( + v + (a — ib)d(\ +2dudv , (2.1) 

where a, b are real constants. Here and throughout this paper a bar will 
denote complex conjugation. The history of the wave is the null hyperplane 
v = and v + = v8(v) where 9{v) is the Heaviside step function, u is a 
second null coordinate. The space-time to the future (v > 0) of the history 
of the wave is Minkowskian and so is the space-time to the past (v < 0) of 
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v — 0. Writing v^2 C = x + iy we see that (271) is in the Rosen-Szekeres form 



ds 2 = -e~ u (e v cosh W dx 2 - 2 sinh W dx dy + e~ v cosh W dy 2 ) +2 e" M dw cfe 

(2.2) 



with 



-u 



-V 



(a 2 + b 2 )v 



+ i 



sinh W 

M 



1 + (a 2 + 6 2 )v 2 - 2av+ 
1 + (a 2 + 6 2 )t; 2 + 2aw + 



1 - (a 2 + 6> 2 ' 
0. 



(2.3) 
(2.4) 

(2.5) 
(2.6) 



We consider the head-on collision of this wave with a wave of similar type. 
This latter wave is described by a space-time with line-element ( |2.2j ) but 
with 



sinh W 

M 



1 - (a 2 + I3 2 )u\ , 

"1 + (a 2 + (3 2 )u\ - 2au. 
1 + (a 2 + p 2 )u 2 + + 2au+_ 
2pu+ 



(a 2 + (3 2 )u 2 + 



2 ' 



0. 



(2.7) 
(2.8) 

(2.9) 
(2.10) 



with a, j3 real constants and u + = u9(u). The history of the wave front is the 
null hyperplane u = in this case. For the collision we consider the space- 
time to have line-element ( |2.2j ) with U, V, W, M given by ( |2.3p - ( |2.6D in the 
region u < 0,v > and given by (|277j)- ( |2.10|) in the region v < 0,u > 0. 
The region u < 0, v < has line-element ( |2.2|) with U = V = W = M = 
(which agrees with ( |2.3| )-( |27l0"|) when both v < and u < 0). The line- 



element in the region u > 0, v > (after the collision) has the form ( f2~2 
with U, V, W, M functions of (w, v) satisfying the O'Brien-Synge || junction 
conditions: If u — 0, v > then [/, V, W, M are given by (|27J)-(|27J) with 
v + = v and if n = 0, u > then U, V, W, M are given by (p77|)— ( pTTOD with 
u + = u. Einstein's vacuum field equations have to be solved for U, V, W, M 
in the interaction region (u > 0, v > 0) after the collision subject to these 
boundary (junction) conditions. These equations are || (with subscripts 
denoting partial derivatives): 



U„ 



(2-11) 
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2V 


= U u V v + U v V u -2(V u W v + V v W u )tanhW , 


(2.12) 


2W UV 


= U U W V + U v W u + 2V U V v sinh W cosh W , 


(2.13) 


2U U M U 


= -2U UU + U 2 + W 2 + V 2 cosh 2 W , 


(2.14) 


2U V M V 


= -2U m + U 2 + W 2 + V 2 cosh 2 W , 


(2.15) 


2M UV 


= -U uv + W u W v + V u V v cosh 2 W . 


(2.16) 



The first of these equations can immediately be solved in conjuction with 
the boundary conditions to be satisfied by U onw = and on v = to yield, 
in u > 0, v > 0, 



(a 2 + b 2 )v 2 -(a 2 + (3 2 )u 2 . 



(2.17) 



The problem is to solve (|2.12|) - ([2.13|) for V, W subject to the boundary con- 
ditions and then to solve ( |2.14j ) and (2.15) for M. Equation ( J2.16|) is the 
integrability condition for (|2.14 ) and ( 2.15|) . 



3 The Backscattered Radiation Fields 



We shall for the moment focus attention on the two field equations ( |2.12|) and 
(|2.13| ). All of our considerations from now on will apply to the interaction 
region of space-time u > 0, v > after the collision. Introducing the complex 
variables 



A = - V u cosh W + iW u 



B = -V v coshW + iW v 



(3.1) 



we can rewrite the two real equations ( p.!2| ) and ( |2.13| ) as the single complex 
equation 

2A V = U u B + U v A - 2iA V v sinh W , (3.2) 
or equivalently as the single complex equation 



2B U = U U B + U v A-2iB V u sinh W . 



(3.3) 

Given the form of the line-element ( |2.2|) it is convenient to introduce a null 
tetrad {m, rh, I, n} in the region u > 0, v > defined by 



m 



e U/2 

~V2 



e v ^ 2 ( cosh — — i sinh — ^ + e v ^ 2 ( sinh — — % cosh — ^ ^~ 
\ 2 2 J dx V 2 2 J dy 

(3.4) 
(3.5) 



n 



d_ 

dv 
d_ 
du 



(3.6) 
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with m the complex conjugate of m. The integral curves of the vector fields 
I and n are twist-free, null geodesies. The coordinate v is not an affine 
parameter along the integral curves of I and these curves have complex shear 
a i and real expansion p\ (we use the standard definitions for these quantities 
given in §4.5 of [ll| for example) given by 

a^-e^S, Pl = \e M l*U v , (3.7) 

with B as in (|3.1|) . Likewise the coordinate u is not an affine parameter 
along the integral curves of n and these curves have complex shear a n and 
real expansion p n given by 

a n = \e M ' 2 A, Pn = 1 -e M ^U Ul (3.8) 

with A as in flO]). We thus see from Q2T71) , ( [OP and (^) that the ratio 



pi/ p n is now known in the region u > 0,v > 0. In terms of the variables 
introduced above the non-identically vanishing scale-invariant components 
|§ of the Riemann tensor in Newman-Penrose notation are \l/o, ^2, ^4 given 
by 

2^ = B V + (M V -U V ) B+iBV v sinhW , (3.9) 
2^ 2 = M uv -- (AB-Ab) , (3.10) 
2^4 = A u + (M u - U u ) A + iA V u sinh W . (3.11) 



When these are non-zero we interpret \l/o as describing radiation, having 
propagation direction n in space-time, backscattered from the wave with 
history u — 0,v > and we interpret ^4 as describing radiation, having 
propagation direction / in space-time, backscattered from the wave with 
history v = 0, u > 0. Thus the integral curves of the null vector fields n 
and / are the 'rays' associated with the backscattered radiation from the two 
separating waves after collision. 

We now look for an interesting assumption to make regarding the rays 
associated with the backscattered radiation fields. Since p\ / p n is known from 
(|3.7p , (|3.8p and (|2.17| ) we focus attention on the complex shears o\ and a n in 
O) and (13.81). Let us write 



A=\A\e ie , B=\B\e^, f = 6 - <f> , (3.12) 
with 9 and real. From (|3.2j ) and (|3.3| ) we can obtain the equations 



\B\ 

o v = — -Aj u u sin f-Vv sinh w > ( 3 - 13 ) 



2\A 
\A 



2\B\ 



U v sin/- V u sinh W , (3.14) 
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and 



2/^ + ^,45-^5) lj sin A - ([/ „ N s i n /) . (3.15) 



The arguments of A, B respectively are tetrad dependent. If we transform 
the tetrad {m, rh, I, n} by the rotation 

m — ► rh = e 1 ^ m , (3.16) 

with ■?/> a real-valued function, then 

0-0 = 0_2V, 0^0 = 0- 2V, /-/. (3.17) 



On account of the field equation Q3.15 ) we can, without loss of generality, 



arrange to have 9 + <fi = constant. To see this we deduce from (|3.13j ) and 
fl3jg ) that 



|4| CA, sin/ + / u -2KsinhW-4^ u , (3.18) 

(^ + ^) = -TTT^ sin/-/„-2KsinhW-4^ • (3.19) 

v /v \A\ 



We are free to choose ip to make the right hand sides of (|3.18 ) and (|3.19|) 



vanish because the integrability condition for the resulting pair of first order 
partial differential equations for ip is the field equation (|3.15 ). Hence it is 



always possible to choose a tetrad so that 9 and in ( |3.12| ) have the property 
that 9 + = constant. This result suggests that to discover an interesting 
assumption to make about the rays associated with the backscattered radia- 
tion fields we should consider the ratio (because it depends upon 9 — 4> and 
not 9 + 4>) 



with the last equality following from (ft.7|) and ( |3.8| ). We note that / satisfies 



the second order equation (|3.15|) to which we shall return later. It is clear 



from ( |3.1|) that a change of parameters u — >• u = u{u) and v — > v = v(v) 
along the integral curves of n and / rescales A and 5 by a function of u and 
a function of v respectively. This change of parameter obviously leaves the 
form of the line-element Q2.2p invariant. Also from the field equation ( |3.2| ) 
we deduce that 

(\A\ 2 ) v -U v \A\ 2 = )-U u (aB + AB) , (3.21) 
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and from the equivalent equation (|3.3| ) we find 

(\B\ 2 ) u -U u \B\ 2 = ]-U V (AB + A 



B 



^From these two equations we obtain 



log 



\A\ 
YB\ 



U„ 



(ab + ab" 
\aT ' 



u„ 



AB + AB 
\B\ 2 



(3.22) 



(3.23) 



which is a partner for the equation ( [3 .15 ) for /. This suggests that it might 
be interesting to explore the following assumption concerning the rays asso- 
ciated with the backscattered radiation: there exist parameters u, v along the 
integral curves of n and I respectively such that \A\ 2 = \B\ 2 (or equivalently 
\a n \ 2 = \o~i\ 2 ). This is equivalent to the assumption that there exist functions 
C(u),D(v) such that 



\A 



\B 



When u — it follows from 



\B\ 



--C(u)D(v) . 

I ) and (|J) that 
4 (a 2 + b 2 ) 



and when v = it follows from 



\A\ 



- (a 2 + b 2 ) v 2 ) 2 ' 
I) and (H) that 

A{a 2 + (3 2 ) 

- (a 2 + f3 2 )u 2 ) 2 ' 



(3.24) 



(3.25) 



(3.26) 



Also when u = we see from ( |2.17| ) that the right hand side of (|3.21|) vanishes 
and thus solving (|3.21| ) for \A\ when u = we obtain 



\A\ 



4(a 2 + /3 2 ) 
1 - (a 2 + b 2 ) v 2 



(3.27) 



with the constant numerator here (the constant of integration) chosen so that 
the two expressions ( |3.26p and fl3.27p for \A\ 2 agree when u = and v = 0. 
Similarly when v — the right hand side of ( |3.22j ) vanishes and we readily 
obtain, when v = 0, 

\ B \ 2 = 1 7 2 T Q2\ 2 ■ ^ 3 ' 28 ) 

1 — (a 2 + p 2 ) u 2 

Thus ( |3.24 ) together with the boundary conditions at u 
results in 



and at v = 



\A[ 
\B\- 



a 



b 2 



{a 2 + b 2 ) v 2 
{a 2 + (5 2 ) u 2 



(3.29) 



Hence there exist parameters (u, v) given by 

sin- 1 (u \ja 2 + f3 2 ^j , v = sin" 1 (y Va 2 + b 2 ) , (3.30) 



u = sin 



such that 

V 2 cosh 2 W + W 2 

V£ cosh 2 W + W 2 ~ 1 ' ' 

We note that when u = we have w = and when v = we have t> = 0. 
We shall express ( |3.31| ) by saying that, in the coordinates (u,v), ( |3.29|) reads 
\A\ 2 = \B\ 2 . In the coordinates (u, v) the form of the field equations and the 
expressions for the Riemann tensor remain invariant, with the derivatives 
with respect to (u, v) being replaced by derivatives with respect to (u, v) and 
with M replaced by M according to 

^/a 2 + P 2 Va 2 + b 2 ' 1 ' 1 

We note from (|2.17|) and ( p.30|) that in the barred coordinates 

e~ u = cos(m — v) cos(m + v) . (3.33) 

Also in these coordinates ( p.23|) becomes 

{Uu cos /) - = (U v cos f). U a fu = U v f v , (3.34) 

(the equivalence here following from ( |3.33| ) since now Una = U m ) from which 
we conclude that 

/ = /(A) with A = CQS ^ " g| . (3.35) 

cos(u + V) 

We see that A = 1 when u = and/or when v — 0. Also ( fj.21 ) and ( |3.22| ) 
become 



J V 


\A\ 2 


= u a 


\A\ 


) -Uu 

1 u 


\A\ 2 


= u v 


\A\ 



(3.36) 
(3.37) 



It thus follows that, again in the barred coordinates, 

\A\ 2 = \B\ 2 = e u g{\) , (3.38) 
for some function g(X) satisfying 

A g' = g cos/, (3.39) 
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with the prime here and henceforth denoting differentiation with respect to 
A. 

In order to intepret physically the implications of our assumption that 
there exists (u,v) such that |v4| 2 = |£>| 2 we proceed as follows: using (|3.20|) 
and the field equations ( |2.14|) and ( |2.15|) in the expressions (|3.9|) and (|3.11|) 



for \&o and ^4 we find that we can write 




-iU v f v -U m + - \B\ 2 + — U U U V + U V (logr^r 



(3.40) 



and 



^4 1 f . TT 1 1 * 1 2 A „ „ „ A \A\ 



U u ^ = -^U u f u -U uu + -\A\ z + — U u U v + U u [\og^J j . (3.41) 

When these are expressed in the coordinates (u, v) we can put \Ay = \B\* 
and as above U m = Um and U a f u = Uy / c and thus 

U *J[ = U ^- ( 3 - 42 ) 



^From fl3.42| ) it follows, using the second of ( |3.7| ) and of ( |3.8|) that 

Pi 2 l*4| 2 = p; 2 |^o| 2 • (3.43) 



This equation, which is a consequence of our basic assumption, deserves some 
physical interpretation. |\l/o| 2 and |^4| 2 are analogous to the energy densities 
of electromagnetic waves propagating in the n and I directions respectively in 



space-time, in the same sense that the Bel- Robinson tensor JTTJ] is analogous 
to the electromagnetic energy-momentum tensor. However with the coordi- 
nates carrying the dimensions of length these quantities have the dimensions 
of (length) -4 . The quantities pj 2 \^\\ and p~ 2 \^o\ 2 both have dimensions 
(length) -2 of energy density and are positive definite expressions in terms of 
the backscattered radiation fields [we note that the backscattered radiation 
has non-vanishing shear and expansion and thus does not consist of systems 
of plane waves]. Hence it seems reasonable to suggest the following interpre- 
tation for the equation ( |3.43| ): in the coordinate system (the barred system) 
in which \A\ 2 = \B\ Z the energy density of the backscattered radiation from 
each of the separating waves after the collision is the same. We note that 
(|3.43j ) also holds for our solution || which describes the collision of an im- 
pulsive gravitational wave with an impulsive gravitational wave sharing its 
wave front with an electromagnetic shock wave. This latter solution contains 
the Khan and Penrose ]I| solution and a solution of Griffiths [[EJ as special 



10 



cases. Also ( |3.43| ) holds (trivially) for the Bell and Szekeres [|T3| solution 
describing the collision of two electromagnetic shock waves. These exam- 
ples demonstrate that ( |3.43| ) holds for a class of collision problems involving 
gravitational impulse waves and/or electromagnetic shock waves. 



4 Integration of the Field Equations 



We begin by writing (|3.15|) in the barred system (u, v). Using ( |3.20| ) with 
\A\ = \B\, ( CT) and ( ggg ) we find that 



;i - A 2 ) f" - 2 A/' + f sin/ = sin/ 



:i-a s 



/'cos/, (4.1) 



with g given in terms of / by ( |3.39| ). We can simplify ( [4. 1| ) to read 



A d 



sin / dX 



;i-a 2 ) r+ 



sin/' 



(4.2) 



We get a single third order equation for / by eliminating g (taken to be non- 
zero) between ( |3.39| ) and (O). Since we are working in the barred coordinate 
system (|3.20|) gives 



A -V n cosh W + iWu 



B 



-Vb cosh W + iWv 



1 - ih 
1 + ih 



(4.3) 



where, for convenience, we have introduced h(X) by the final equality. After 
eliminating g from (|3.39|) and ( |4.2|) it is useful to write the resulting equation 
differential equation for h(X). Then defining 



G 



X 2 ) h! - 



li 



1 + h 2 

the equation for h(X) can be put in the form 

4 



(4.4) 



XG" + G' - - G 
X 



GQ 
1-A 2 



where 



Q = ^(l-h 2 )G' + ±-(l + h 2 )G . 



We remark that if we define 

A 



1 



P=-(l + h*)G' + -(l-h 2 )G 



(4.5) 
(4.6) 

(4.7) 
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then 



XP' = Q 



and 



P 2 -Q 2 = \ 2 (Gf-AG 2 



(4- 



In studying the differential equation (|4.5| ) for h we found it helpful to write 
and the second of (14.80 in the form 



A G" + G' 



4 
A 

P 2 - A 2 (P 



G 



/\2 



XGP' 

T^x 2 ' 

V 2^2 . 4G 2 



A (G' 



(4.9) 
(4.10) 



and to work with these equations. Before proceeding further however we 
need to know h and h! when u = and/or when v — 0, i.e. we require /i(l) 
and fo'(l). To find start by writing (|4.3|) , using ( |3.30|) , as 



l-ih{\) _ Va 2 + b 2 ^Jl-(a 2 + P 2 )u 2 / -V uCOS hW + iW u 
l + ih(X) ~ ^/ a 2 + ( a 2 + &2 ) v 2 \-V v coshW + iW V/ 



(4.11) 



and evaluate this equation when u = and v = 0. From the boundary 
conditions on V and W given by (|2.4|) , (|2.5|) , (|2~8|) and (|2T9| ) we have that 
when u = and n = 0: 



K = 2« ,K = 2a ,W U 



2 /?, W„ = 2 b 



and thus from ( [4.11] ), 



where 



1 -ih(l) _ y/a 2 + b 2 fa-i/3" 
l + ih(l) ~ ^a 2 + f3 2 \a-ib / 

ia <X-iP 



V« 2 + P' 2 

It thus follows from (f4.13|) that 



Mi) 



tan 



'a — p 



a — ib 
y/a 2 + b 2 ' 



k (say) 



(4.12) 



(4.13) 



(4.14) 



(4.15) 



Next to find h'(l) we begin with (|4.11|) and by two differentiations obtain 
from it 



U {a 2 + (5 2 )h'{\) 

(i + zMi)) 2 



d 2 ( -V u coshW + iW u 



dudv \ —V v cosh W + iW v 



(4.16) 



. (u=0,v=0) 



To evaluate the right hand side here we first note from (|2.17[) that when 
u = and v = 0, U u = and U v = 0. Also from the boundary conditions on 
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W we have W = when u = andv = 0. Now evaluating the field equations 



2.12j) and (|2.13|) when k = and i> = we easily see that in this case 

V,. 











(4.17) 



^From the boundary conditions satisfied by V and W we have, when u = 
and v = 0: 

V w = V uu = W vv = W uu = . (4.18) 

Next differentiating ( |2.12|) and ( 2.13|) with respect to u we find that when 
u = and v = 0: 



2 a 2 a — 6 /3 2 a — 8 b a (3 , 



2b(a 2 + f3 2 ) + 8aaf3 . (4.19) 



Finally differentiating ( p,12j ) and (p.!3| ) with respect to v we find that when 
u = and v = 0: 



V' 



2a a 2 - 6a 6 2 -Sab (3 



2f3(a 2 + b 2 ) + 8aab . (4.20) 



Now substituting all of these results into the right hand side of ( |4.16|) we 
obtain 

ti(l) = k , (4.21) 
with k given by ( |4.15| ). Using ( |4.15| ) and ( |4.21|) in (^4|) we see that 

G{1) = = • (4.22) 

We can now set about solving ( |4.9| ) and ( |4. lQ ) for G and then obtain h(X) 
from (ELI]) . 

Differentiating ( |4 . 1 0| ) with respect to A and using ( |4.9| ) we find that either 
(a) p> = o or (b) if P' ^ then 



x p// , pi 1 p _ AGG/ 



(4.23) 



We can quickly dispose of case (a). If P = constant ^ then ( |4.10|) can be 
integrated to yield 

1 \ 



2 \±2 



c§A±* 



(4.24) 



where Co is a constant of integration. It is easy to see that this constant 
cannot be chosen to satisfy both boundary conditions (4.22). Also if P = 
then (|4. 10|) integrates to 

G = Cl A ±2 , (4.25) 
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with ci a constant of integration. Clearly we must have c± = to satisfy 



(|4.22j ). Thus the only acceptable solution in case (a) is G = 0. Turning now 
to case (b) with fl4.23| ) holding we find that we can integrate this equation 
once (using ( |4.22|) ) to read [fl4|| 




AG 2 



2(1 -A 2 ) 

^From the first of (|4.8|) this can be written 



(4.26) 



(P + Q) + A 2 (P-Q) = ^G 2 , (4.27) 



and from (14.71) and (14.81) this reads 



Ar / Afc 2 (1 - A 2 fr 2 ) 

A glance at ([4.3|) and ( |4.4| ) shows how h and thence G are constructed from 
the functions V and W appearing in the line-element (|2.2|) for u > 0, v > 
0. On the boundaries of this region A = 1 and within this region A > 1 



with A becoming infinite when the right hand side of ( 2.17 ) vanishes. The 
A = constant > 1 curves densely fill the interior of the region «S(say) with 
boundaries (&i) u = 0,v > 0, (6 2 ) v = 0,u > and (6 3 ) the right hand 
side of ([2.171 ) vanishing with u > 0,v > 0. Within S there is one curve 
A = constant > 1 passing through each point. When the field equations 
are completely integrated for (u, v)eS the boundary (6 3 ) turns out to be a 
curvature singularity. For G analytic in S we conclude from (|4.22|) and ( |4.28| ) 



that G = in S. It thus follows from (|4.4| ) with the boundary condition (|4.15|) 
that 

h(\) = k\, (4.29) 

for (u, v)eS. 

We are now at the following stage in the integration of the field equations: 
the function U in the line-element fl2.2|) is given by ( |2.17|) in coordinates (u, v) 
or by ( |3.33|) in coordinates (u,v). Also on account of ( |4.3|) and (|4.29|) the 



functions V and W in ( |2.2| ) satisfy the differential equation 

A -14 cosh W + iW a _ 1 - i k A 
B ~ -Vv cosh W + iW„ ~ 1 + ikX ' 



(4.30) 



with A given by ( p.35|) and k by ( |4.15| ). We shall now solve this complex 



equation for V and W in terms of the barred coordinates. First we need to 
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note the boundary values of V and W in terms of the barred coordinates. 
By (|2.4j) and ( p.30| ) we have when u — 0, 



-\/a 2 + b 2 — a sin v ) +6 2 sin 2 v 



V^ 2 + 6 2 + a sin v ) + fe 2 sin 2 t> 



1 — e J/3 sin v 



1 + e^ 3 sin v 



(4.31) 



with the second equality following from ( f4.14| ). By ( |2.5| ) and ( |3.30| ) we have 
when u = 



sinh 



2 6 sin f 



(e ^ sin f — e^ sin f 



-y/a 2 + b 2 cos 2 w 



smv 



(4.32) 



The corresponding boundary values on v = are obtained by replacing u by u 
and /3 by a in the final expressions in (|4.31| ) and CN= . 3 2|) . It is convenient to use 
a complex function E (the Ernst function) in place of the two real functions 
V and W defined (in a way that is suggested by the final expressions in (f4.31|) 
and (|Q2D ) by 



-v 



(1-E) [1-E 



sinh W 



E-E 



1 - \E 



' _(1+E) (l + E) 
or equivalently by 

sinh V cosh W + i sinh W 

E = 

1 + cosh V cosh W 

Now (|4.31|) and ( |4.32|) can be written neatly as: 
when u = , E 

and correspondingly 

when v = , E 

In terms of £7, the complex functions A, B can be written 

2 cosh W 



2 ' 



(4.33) 



— 7/3 • - 

e H sin v 



e * a sin u 



A 



2 cosh W - 
E,-, 



B 



1-E 2 ' 1-E 2 
Substitution into ( |4.3(J| ) simplifies this equation to 



Er, 



(4.34) 



(4.35) 



(4.36) 



(4.37) 



En, 



Eu = ik\ (E n + Eg) . 



(4.38) 
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With A given by (|3.35| ) this equation establishes that 

E = E(w) with w = sm(u + v) + ik sinfu — v) . 



(4.39) 



We can now determine E using the boundary conditions ( 4.35Q and ( [4.36 ). 
To see this easily we first write k in fl4.15| ) in the form 



( 



-if) 



(4.40) 



Using this in ( [4.39 ) we see that we can consider E to have the functional 
dependence: 



E = E (e~ ta sin u cos v + e~ l/3 cos u sin v 
Now the boundary conditions ( |4.35 ) and ( i.'St ) establish that 
E = e~ ia sin u cos v + e _l/3 cos u sin v , 



(4.41) 



(4.42) 



and thus the functions V, W appearing in the line-element (2.2) are deter- 
mined by ( |4.33| ) in the coordinates (u, v). They are then converted into the 
coordinates (u,v) using the transformations ( p.30|) (see section 5 below). 

Finally in the barred system the field equations (|2.14j) and (|2.15f) for M 
read (using ( |4.33j ) and ( |4.37|) ) 



Mr, 



u, 

~U r , 



2 E n E a 



Uu (l - |£| 2 ) 



2 > 



1 2 EyEy 

2 U v (l-\E\ 



2 > 



(4.43) 
(4.44) 



with M related to M by ( |3.32j ). Since we must have M = when u = and 
when f = we see from ( |3.32| ) that 



and 



when u = 



when u = 



-JW 



-M 



COSf 



cosu 



(4.45) 



(4.46) 



In ( |Q3D and JO^) , £7 is given by 
find 



Eu E u 



\E\ 



V« 2 + f3 2 Va 2 + b 2 ' 
) and E by flO^ ). Using ( ^ ) we 



(4.47) 
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with 



1 - \E 



2 2 
= COS 



a — p 



cos 2 (u + v ) + sin 2 



a — p 



cos 2 (w-v). (4.48) 



The only complication in solving (|4.43 ) and ( 4.44 ) is in dealing with the final 
term in each. In the case of (|4.43D this now involves evaluating the integral 



2du 



Uu 1- |J5?| 



2\d\ 



(A 2 - 1) {cos 2 + A 2 sin 2 



(4.49) 



where we have changed the variable of integration from u to A, given in 
( |3.35| ), with v held fixed. This integral is easy to evaluate and using Q4.48|) 
again we obtain from ([4.43 ) 



1 - \E\ 



F(v) Jcos(u — v) cos(m + v) 



(4.50) 



with F(v) a function of integration. By ( (4.45[) we find that in fact F is a 
constant given by 

(4.51) 



a 



2 + (3 2 Va 2 + b 2 . 



It is straightforward to see that ( 4.50|) with ( |4.51| ) also satisfies ( |4.44| ). The 
integration of Einstein's vacuum field equations is now complete. 



5 Discussion 

The purpose of this paper has been to propose a simple key to open up the 
boundary-value problem which is involved in deriving a model in General 
Relativity of the vacuum gravitational field left behind after the head-on 
collision of two plane impulsive gravitational waves. This 'key' has been 
provided, with some motivation, in section 3 (following equation (|3.23| )) and 



a physical interpretation has been suggested there for the interesting equation 
(|3.43| ), which is a consequence of the key assumption and some of the vacuum 
field equations. Our approach has been to focus attention on properties of 
the backscattered gravitational radiation present after the collision. This 
appears to be a non-linear phenomenon whose presence therefore ought to 
be expected to play a central role in the development of a scattering theory 
for gravitational radiation. 

Notwithstanding the simplicity of our key assumption, the derivation of 
the line-element of the vacuum space—time in the region S in section 4 [the 
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region S is defined following equation ( [4.28|) 1 is complicated. It is therefore 



useful to summarise the result: the vacuum space-time in the region S has 
line- element of the form Q2.2p with U given by ( |2.17| ), V and W given by 
(|4.33| ) with the complex function E in ()4.42|) expressed in coordinates (u, v) 



as 



E=(a + iP)uyJl - (a 2 + b 2 ) v 2 + (a + ib) v yjl - (a 2 + (3 2 ) u 2 , (5.1) 



while M is reconstructed in coordinates (u,v) using (|3.30|) , ( |3.32| ), (|4.50|) and 
( |4.51| ). The result is 



(a 2 + (3 2 ) u 2 yjl - (a 2 + b 2 ) v 2 ' 



with e u given in ( PT7| ). If in ( gT7|) , ( |5TT| ) and (^) we put a 2 + b 



a: 



+ (3 2 we recover the original form of the Nutku and Halil solution. 
If in addition b = (3 = (and thus E = E and so W = 0) we recover 
the original form of the Khan and Penrose QXJ solution. We note that in the 
region S a curvature singularity [|TJ , H is encountered on the boundary where 
(a 2 + b 2 ) v 2 + (a 2 + (3 2 ) u 2 = 1 and the solution above is valid only up to this 
space-like subspace. 

Finally we wish to emphasise that the approach to solving collision prob- 
lems in General Relativity developed in this paper is not confined to the ex- 
amples worked through above. With different boundary conditions to those 
described by equations ( |2.3|) -( [2~T0|) and the paragraph following ( |2.10| ), the 



technique is capable of solving new collision problems (see, for example ||) 
This collaboration was funded by the Department of Education and Sci- 
ence and by the Ministere des Affaires Etrangeres. 
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This follows because ( |4.23|) can be written ^ Q — A J X 



where 



X = A P'+ (fz^s) -P— 2H • The constant of integration then vanishes 
on account of (4.22), which in particular implies by ( (4.7[ ) that -P(l) = 0. 



A Plane Impulsive Wave 

The line-element of the vacuum space-time describing the gravitational 
field of an impulsive pp-wave || is given by 

ds 2 = -2dZ dZ + 2dV (dU + 5(V) f (z, 2) dv) . (A.l) 

Here / is a real- valued function which is harmonic in Z, Z: 

d 2 f 

-^4- = • (A.2) 
OZdZ K 1 
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The only non-vanishing component of the Riemann tensor in Newman- 
Penrose notation is 

%=^S(V), (A.3) 

which is therefore Petrov type N with d/dU as degenerate principal null 
direction. Integral curves of d/dU, which have vanishing expansion and 
shear, generate the history of the wave-front V = 0. Hence V = is a 
null hyperplane. Following the example of Penrose J7| a discontinuous coor- 
dinate transformation removes the 5- function from the line-element (|A.l ) 



and introduces a coordinate system (v,u, in which the metric tensor is 
continuous (C°) across V = 0. This transformation is given by 

V = v, U = u-6(v)f{(,0 + \g\ 2 v + , Z = C + v + g(C), (A.4) 

where 6{v) is the Heaviside step function, v + = v9(v) and g(() = —df/d( is 
an analytic function of (. The line-element ( |A.1|) is transformed under (|A.4|) 
into the Rosen form 



ds 2 



d( + v + K(()d( +2dudv, (A.5) 



with K(C) = dg/d( = -d 2 f/dQ 2 and O becomes 

9o = j^6(v). (A.6) 

For a plane impulsive wave with two degrees of freedom of polarisation 
d 2 f/d( 2 is a complex constant. We take / = — Re {(a + ib) ( 2 } with a, b 
real constants for this case. For a linearly polarised plane impulsive gravita- 
tional wave either a = or b = 0. We note that ( A.4 ) incorporates Penrose's 



geometrical construction of the impulsive pp-wave whereby the history of the 
wave is formed in Minkowskian space-time by first subdividing the space- 
time into two halves v > and v < each with boundary v = and then 
reattaching the halves on v = identifying the points (v = 0, u, () and 
(v — 0, u — /(£, C))0- This is a mapping of the two copies of v — where 
points on the same generators £ = constant of v — are mapped one to 
the other by the translation u — > u — /(C, C) an d the mapping preserves the 
intrinsic degenerate metric on v — 0. 



B Incoming Waves Linearly Polarised 



If the incoming waves are linearly polarised b = in (|2.3|) - (|2.5|) and /3 = 
in (|2.7|)-(PTP|). In the interaction region (u > 0, v > 0) after the collision the 
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function W = 0. This is the problem solved by Khan and Penrose Our 
basic assumption following ( |3.23|) in this case reads: there exist parameters 
(u,v) along the integral curves of n and I respectively such that V u = Vy. 
These parameters are given by ( |3.30| ) with b — (3 — 0. Now V = V(u + v) in 
u > 0, v > 0. In the barred coordinates the boundary value of V when u = 
is obtained from (|2.4[ ) with b = to be: 

e v = 1 + Sin " . (B.l) 
1 — sin v 

Thus for u > 0, v > we have 

v 1 + sin(w + v) (cosm + sin w) (cos w + sin-u) 

e = . , = t z ^—77 z — t ■ (B.2) 

1 — sm(u + (cos u — sin dj (cos v — sin tt) 



Using (|3.30| ) with b = (3 = we can write ( p.2| ) as 



_ (vl - oi 2 u 2 + aw)(Vl - a?v 2 + an) 
(yl-aV — au)(\/l — a 2 f 2 — aw) 

which is the Khan and Penrose [Q expression for V in the interaction region. 
From (|3.30|) and ( p.33|) we have U in the coordinates (u, v) given by 

e~ u = 1 - a 2 w 2 - a 2 t» 2 . (B.4) 

The only remaining function to be determined is M. In the coordinates 
(u, v) this is replaced by M with the latter related to M via (|3.32| ) with 
b — (3 — 0. If we let Q — M + then in the barred coordinates the field 
equations ( [2.14| ) and ( |2. 15| ) reduce to the following equations for Q: 

Q a = Q- = 2 tan(u + v) . (B.5) 

Since M = when u = or v = it is easy to see that when u — 0, 
Q = log (cos 2, u/aa) and when -y = 0, Q = log (cos 2 w/a a). It thus follows 
from (|B~5D that 

0= /cos^n 



Hence by ([H|) with 6 = /3 




a a 



—M __aa_ -Q+iu [cos(u + v) cosju - v)}? 

e = z ~e * 2 = z z Tfi zr- • (B.7J 

COS U COS V COS U COS V cos^ [u — v) 

In the (m, f ) coordinates this reads 

e-» = , ( 1 -^ 2 - aV ) f , fB .8) 



Vl — a 2 u 2 Vl — a 2 f 2 [-\/l — ot 2 u 2 yl — a 2 ^ 2 + aauv] 
which is the Khan and Penrose M expression. 



2 



21 



